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1.  INTRODUCTION 


Over  the  past  few  decades  microwave  Doppler  radars  have  demonstrated 
the  ability  to  measure  wind  fields  in  regions  where  there  is  sufficient 
reflectivity  from  particles,  such  as  precipitation,  to  provide  a  radar  echo. 
Recent  experiments  have  shown  that  sensitive  Doppler  radars  can  detect 
and  measure  the  radial  velocity  of  refractive-index  fluctuations  in  optically 
clear  air.  Using  the  Velocity-Azimuth  Display  (VAD)  technique.  Frequency 
Modulated-Continuous  Wave  (FM-CW)  Doppler  radar  and  pulse  Doppler  radar 
can  measure  average  wind  profiles  from  either  precipitation  scatterers  or 
refractive- index  scattering.  Therefore,  microwave  Doppler  radar  is  capable 
of  providing  wind  measurements  in  all  weather  conditions.  However,  a 
vertical  profile  of  the  average  horizontal  wind  near  a  particular  point, 
such  as  the  touchdown  point  or  middle  marker,  may  not  be  adequate  to  warn 
aircraft  of  the  presence  of  dangerous  winds  during  take-off  and  landing 
when  the  wind  field  Is  spatially  or  temporally  variable.  The  radar  must 
be  capable  of  accurately  measuring  the  radial  velocity  at  low  elevation 
angles  and  short  ranges  to  determine  the  wind  profile  encountered  by 
aircraft  on  the  landing  and  take-off  paths  near  the  runway.  An  experi¬ 
mental  program  is  being  conducted  at  the  National  Aviation  Facilities 
Experimental  Center  (NAFEC)  to  determine  if  an  air  traffic  control 
radar,  equipped  with  a  suitable  antenna  and  with  Doppler  data  processing, 
has  sufficient  sensitivity  and  resolution  to  accurately  measure  radial 
wind  profiles  at  low  elevation  angles  and  short  ranges  In  optically 
clear  air  (as  well  as  In  regions  of  precipitation). 

The  mean  Doppler  velocity  depends  on  the  radial  velocity  distribution 
of  the  scatterers  In  the  radar  resolution  cell,  so  that  the  radar -measured 
velocity  Is  a  filtered  version  of  the  actual  wind  component  along  the 
axis  of  the  radar  beam.  Velocity  gradients  existing  over  distances  much 
smaller  than  the  radar  beamwldth  will  cause  a  broadening  of  the  Doppler 
spectrum  but  will  not  be  revealed  by  the  radar-measured  mean  velocity. 
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Velocity  gradients  exlstlnq  over  distances  comparable  with  the  beaimvldth 
will  be  smoothed  so  that  peak  values  will  not  be  detected.  Furthermore, 
If  the  radar  Is  offset  from  tlie  runway  so  that  the  axis  of  the  radar 
beam  Is  parallel  to  but  displaced  from  the  flight  path,  the  measured 
wind  component  will  be  representative  of  the  actual  wind  component  along 
the  flight  path  only  If  the  wind  field  Is  horizontally  homogeneous  (i.e.. 
If  there  are  no  horizontal  gradients). 

This  report  describes  a  computer  simulation  for  testing  the  effects 
of  both  finite  radar  resolution  volume  and  radar  displacement  relative  to 
the  runway  on  the  radar-measured  component  of  the  wind  along  the  flight 
path.  The  simulation  program  contains  six  wind  shear  models,  and  the 
results  obtained  with  these  models  are  Included  In  this  report.  The 
simulation  program  is  general  enough  to  acconmodate  any  analytical  wind 
model  that  expresses  the  Cartesian  wind  components  In  terms  of  Cartesian 
position  coordinates. 


2. 


FORMULATION  OF  THE  SIMULATION  PROGRAM 


The  location  of  the  radar  and  the  runway  are  Illustrated  in  Figure 
2.1.  The  nominal  touchdown  point  is  the  origin  of  the  Cartesian  coordi¬ 
nate  system.  The  radar  location  is  at  (XR,  YR,  0).  The  azimuth  pointing 
angle  (AZ)  is  measured  clockwise  from  the  positive  y  direction,  and  the 
elevation  pointing  angle  (EL)  Is  measured  upward  from  the  horizontal. 

The  aircraft  travels  alonq  the  glide  path  at  an  angle  $G  from  the  horizontal, 
so  that,  at  a  radial  distance  RG  from  touchdown,  its  coordinates  are 
given  by 


X  -  PG  cos  *G 
y  «  0 

z  •  Rg  sin  *G. 

The  analytical  models  of  the  wind  field  provide  the  Cartesian  wind 
components  u,  v,  and  w  for  any  (x,y ,z).  The  wind  component  in  the  direc¬ 
tion  of  the  aircrafts'  flight  path  is  calculated  from 

VG  ■  u  cos  $G  ♦  w  sin  $G. 

With  this  formulation,  a  tailwind  is  a  negative  qlide-slope  velocity  and 
a  headwind  is  a  positive  velocity  for  an  aircraft  landing. 

A  spherical  coordinate  system  with  Its  origin  at  the  radar  has 
coordinates  (r,-,o)  that  are  related  to  the  Cartesian  coordinates  by 

x  ■  r  sine  coso  ♦  XR 
y  *  r  cosO  coso  ♦  VR 
z  ■  r  sino. 

The  radial  component  of  the  wind,  at  any  point.  Is  related  to  the  Cartesian 
wind  components  by 
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Figure  2.1  Radar  location  and  orientation  relative  to  the  glide  path. 


VR(r,0,$)  =  u(r,0,ij>)sin9  cos$ 

♦  v(r,0,$)cosd  cos$ 

+  w(r,e,t)s1n4>. 

We  assume  uniform  radar  reflectivity  and  calculate  the  theoretical  radar- 
measured  radial  velocity  at  the  point  (R,AZ,EL)  from 

/  VR(r.0,*)W(r,e,0)dV 

Vm(R,AZ.EL)  «  -  . 

/  W(r.6,*)dV 

V 

where  W  is  the  two-way  radar  antenna  Illumination  function  defined  by 

w .  „pi  , . 

2a 

g 

and  A  *  ( -  ),  where  B  is  the  1-way  antenna  beamwidth.  The  volume 

element  dV  *  (dr)(rdi)(rcosi>dO).  The  limits  of  the  Integration  are  taken 
as  r-R  <  ct/4,  1 0-AZ ,  <  mA,  and  jv?-EL  <  mA  where  t  is  the  radar  pulse 
duration,  c  is  the  velocity  of  propagation  and  m  is  a  constant  (equal 
to  3.1)  which  assures  that  the  error  introduced  by  iqnoring  part  of  the 
antenna  illumination  is  negligible.  The  integration  extends  beyond  the 
-20  dB  illumination  points.  The  integral  is  evaluated  numerically, 
using  volume  elements  of  the  order  of  10  m  .  The  actual  wind  component 
along  the  glide  path  and  the  theoretical  radar-measured  radial  velocity 
are  displayed  for  comparison.  The  calculated  radial  velocities  are 
expected  values  that  would  be  measured  by  an  actual  radar  with  random 
measurement  error.  The  radar  resolution  volume  low-pass  filters  the 
true  radial  wind  field  so  the  resolution  of  the  measured  wind  depends 
on  the  pulse  length  of  the  radar  and  the  antenna  beamwidth. 
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3.  WIND  SHEAR  MODELS 


Six  wind  shear  models  were  included  in  the  simulation  program.  Models 
1  through  5  do  not  contain  horizontal  gradients  and  are  therefore  vertical 
profiles  of  the  horizontal  wird.  These  profiles  are  the  same  at  every 
(x,y)  location.  The  first  four  profiles  were  developed  by  Flchtl  and 
Camp1  of  NASA’s  Space  Sciences  Laboratory  for  the  FAA.  The  fifth  profile 
was  presented  In  a  memo  by  Langwe11i.  The  sixth  model  was  developed  by 
the  authors.  It  Is  a  three-dimensional  wind  model  with  horizontal  and 
vertical  gradients  of  all  three  wind  components. 

Profile  1  (Flchtl  and  Camp’s  profile  1)  is  a  logarithmic  profile 
that  Is  typical  of  the  wind  in  the  neutral  boundary  layer.  The  functional 
form  Is 


V(Z)  -  V(20)  )  for  1.5  <  Z  <  1500  feet. 


where  V ( Z )  Is  the  magnitude  of  the  horizontal  wind  at  height  Z  and  V(20) 
Is  the  reference  level  wind  taken  as  11.5  m/s.  The  wind  direction  Is 
constant  with  height  and  the  vertical  wind  Is  zero. 


Profile  2  (Flchtl  and  Camp’s  profile  3)  Is  the  recormended  frontal 
wind  profile.  The  profile  Is 


V(Z)  •  3.5(Z/20)°’43(m/s)  0  <  Z  <  1500  feet 


e  - 

0 

6  -  0.4(Z-100)*9o 

e  •  120  ♦  e 

0 


0  <  Z  <  100  feet 
100  <  Z  <  400  feet 
400  <  Z  <  1500  feet. 


*  G.  H.  Flchtl  and  D.  W.  Camp,  letter  to  F.  Melewlcz  dated  January  23,  1976. 
1  L.  Langwell,  FAA  memo  dated  February  17,  1976. 
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i)  is  the  direction  (in  degrees)  toward  which  the  wind  is  blowing,  measured 

counterclockwise  from  the  x-axis,  and  6  is  the  wind  direction  at  the 

o 

surface.  The  vertical  wind  is  zero. 

Profile  3  (Fichtl  and  Camp's  profile  2)  is  recommended  for  representing 
the  nighttime  stable  boundary  layer.  The  wind  direction  is  constant 
with  height.  The  magnitude  of  the  horizontal  wind  is  given  by  a  linear 
interpolation  from  the  following  table. 


Z(feet) 

V( knots) 

30 

2.7 

75 

4.1 

150 

6.9 

300 

8.3 

450 

11.4 

600 

10.5 

750 

9.2 

900 

8.2 

1050 

8.0 

1200 

7.9 

1300 

7.8 

1500 

7.8 

Below  30  feet  the  wind  is  given  by  V  a  2.7  (Z/30)’^.  The  vertical  wind 
is  zero. 

Profile  4  (Fichtl  and  Camp's  profile  4)  represents  cold  air  outflow 
from  a  thunderstorm.  The  profile  shows  the  winds  encountered  by  the  air¬ 
craft  at  various  heights.  It  Is  a  one-dimensional  representation  that 
does  not  include  the  horizontal  gradients  that  would  be  found  with  thunder 
storm  outflows.  The  horizontal  wind  components  are  given  in  the  following 
tible. 
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The  vertical  wind  in  ft/sec  is  given  by 


w  -  (-5.04  )sin(  ”J)  ') 


iO  *  Z  •  330  feet 


nn  ^w460-2wv460-Zv2  .  ,,460-Zu590-7, 

(1Q.13){  )  -H*2?.ll(  1Td-)(  Y5Q-), 


w  *  -16.94  sinr ' 


330  7  <■  460  feet 

460  7  •  770  feet  . 


Profile  5  (langweil),  an  alternative  to  profile  3,  his  been  IhOwn 
to  be  more  difficult  for  pilots  to  handle  in  flight  trainers.  It  1» 
given  by  linear  interpolation  from  the  following  table. 


V {knots ) 
7.7 
11.6 

19.6 

23.7 
32.4 
30 

26.2 

24.6 
22.9 

22.6 
22.3 
22.3 


HfeetJ 

30 

75 

150 

300 

450 

600 

750 

900 

1050 

1200 

1300 

1500 


Wind  direction  is  constant  with  height  and  there  is  no  vertical  wind. 


Wind  shear  model  6  is  an  analytic  model  that  simulates  a  downburst 
from  a  thunderstorm  as  depicted  by  Fujita  and  Caracena’.  The  center  of 
the  downburst  (XC*YC)  (see  Figure  3.1)  is  adjustable.  The  direction  of 
the  horizontal  wind  is  radially  outward  from  the  downburst  center, 
and  the  magnitude  is  given  by 


VD  -  1502  e 


•oozr 


)  m/$. 


where  Z  is  the  altitude  (km)  and  R  is  the  distance  from  the  downburst 
center  (km).  The  horizontal  wind  components  are 


u  •  VR(X-Xc)/R  m/s 
v  *  VR(Y‘Yt)/R  m/S- 


*  T.  T.  Fujita  and  F.  Caracena,  "An  analysis  of  three  weather -re la ted 
aircraft  accidents",  in  dull.  An.  Met.  Soc.,  58,  1977,  pp.  910-918. 
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The  maximum  horizontal  winds  occur  560  m  from  the  downburst  center.  The 
vertical  wind  is  given  by 


In  the  computer  simulations  a  minimum  radar  range  of  450  m  was 
used.  This  minimum  range  was  chosen  so  that  the  location  of  the  radar 
could  be  taken  as  0  m  above  the  runway  without  having  negative  heights 
in  the  pulse  volume.  The  minimum  height  for  calculating  radar  measured 
winds  was  23.5  m.  The  minimum  range  of  actual  radars  Is  expected  to  be 
at  least  450  m  because  of  ground  clutter. 


11 


4.  SIMULATION  Rf SUITS 


Profiles  1-5  have  no  horizontal  position  dependence;  consequently, 
the  offset  position  of  the  racar  relative  to  the  runway  Is  immaterial 
for  wind  speed-height  profile*.  The  wind  direction  for  profiles 
1,  3,  and  5  was  made  constant  with  height  and  parallel  to  the  runway. 

For  profile  2,  the  wind  direction  varied  with  height  and  the  direction 
at  the  ground  was  also  taken  parallel  to  the  runway.  Figure  4.1  through 
4.15  show  the  simulation  results,  with  three  figures  for  each  profile. 

The  actual  wind  profile  encountered  by  the  aircraft  Is  shown  as  a  solid 
line;  the  simulated  radar-measured  points  are  shown  as  diamonds.  The 
first  figure  in  each  group  of  three  Is  a  wind  speed-height  profile  for  a 
3-degree  antenna  elevation  anqle,  and  the  third  fiqure  in  the  group  is 
a  wind  speed-height  profile  for  a  6-deqree  elevation  angle.  Since  the 
offset  position  of  the  radar  Is  immaterial.  It  was  taken  as  zero  (radar 
located  at  touchdown).  Therefore,  the  radar  ranqe  and  the  distance  along 
the  glide  slope  are  the  same  only  for  a  3-degree  elevation  anqle.  The 
second  figure  In  each  group  shows  the  winds  along  the  3-degree  glide  slope. 

In  all  cases,  the  minimum  ranqe  of  the  radar  was  450  m  (3  microseconds), 
which  corresponds  to  a  minimum  altitude  of  23  m  for  a  3-degree  elevation. 

For  profiles  2  through  5,  the  input  wind  profile  extends  to  a  1500-foot 
altitude;  above  1500  feet,  the  computer  program  set  the  wind  to  zero. 

This  causes  an  unrealistic  wind  shear  In  the  input,  but  It  serves  to 
Illustrate  the  problem  of  finite  antenna  heamwidth  at  long  range  as 
discussed  below.  For  profile  1,  the  input  model  does  not  return  to  zero 
at  1500  feet. 

Figures  4.1  through  4.26  present  simulation  results.  The  legend 
appearing  below  contains  the  information  comon  to  every  figure  and  the 
variables  specifically  appropriate  to  each. 
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Legend 


The  following  parameters  are  common  to  Figures  4.1  through  4.26. 

Aircraft  glide  slope  -  3  deg  Minimum  radar  range  -  450  m 

Runway  azimuth  angle  -  90  deg  Radar  pulse  length  -  1  microsecond 

Radar  beanwidth  -  1.5  deg  Sampling  Increments  -  10  m 

All  plotted  wind  speeds  are  in  meters  per  second 

The  following  table  lists  the  variables  specifically  appropriate  to  each  figure. 


FIG.  NO. 

AZ 

EL 

XR 

YR 

Profile 

e 

XC 

YC 

4.1 

90 

3 

0 

0 

1 

90 

- 

- 

4.2 

90 

3 

0 

0 

1 

90 

- 

- 

4.3 

90 

6 

0 

0 

1 

90 

- 

- 

4.4 

90 

3 

0 

0 

2 

90* 

- 

- 

4.5 

90 

3 

0 

0 

2 

90* 

- 

- 

4.6 

90 

6 

0 

0 

2 

90* 

- 

- 

4.7 

90 

3 

0 

0 

3 

90 

- 

- 

4.8 

90 

3 

0 

0 

3 

90 

- 

- 

4.9 

90 

6 

0 

0 

3 

90 

- 

- 

4.10 

90 

3 

0 

0 

4 

90 

- 

- 

4.11 

90 

3 

0 

0 

4 

90 

- 

- 

4.12 

90 

6 

0 

0 

4 

90 

- 

- 

4.13 

90 

3 

0 

0 

5 

90 

- 

- 

4.14 

90 

3 

0 

0 

5 

90 

- 

- 

4.15 

90 

6 

0 

0 

5 

90 

- 

- 

4.16 

90 

3 

0 

0 

6 

- 

.5 

.5 

4.17 

90 

3 

0 

0 

6 

- 

.5 

.5 

4.18 

90 

6 

0 

0 

6 

- 

.5 

.5 

4.19 

90 

3 

-0.1 

-0.1 

6 

- 

.5 

.5 

4.20 

90 

3 

-0.1 

-0.1 

6 

- 

.5 

.5 

4.21 

90 

6 

-0.1 

-0.1 

6 

- 

.5 

.5 

4.22 

90 

3 

-2.5 

0.1 

6 

- 

.5 

0 

4.23 

90 

3 

-2.5 

0.1 

6 

- 

.5 

0 

4.24 

90 

6 

-2.5 

0.1 

6 

- 

.5 

0 

4.25 

60 

3 

0 

0 

2 

90* 

- 

- 

4.26 

120 

3 

0 

0 

2 

90* 

• 

• 

•  Wind  direction  at  the  ground. 
All  angles  in  degrees 
All  distances  in  km 
d  ■  wind  direction 
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Profile  1  (Figures  4.1  through  4.3)  has  very  small  gradients,  so 
the  radar  has  no  problem  measuring  this  profile. 

Profile  2  (Figures  4.4  through  4.6)  illustrates  the  limitations 
imposed  by  antenna  beamwidth.  The  large  gradients  near  the  100-m 
altitude  are  represented  reasonably  well  by  the  radar  measurements,  but 
the  gradient  at  the  top  of  the  input  profile  is  not.  This  is  because  the 
ranqe  at  the  top  of  the  profile  is  approximately  9  km  and  the  beamwidth 
is  about  235  m.  The  gradient  at  1500  feet  is  smoothed  by  the  235-m 
vertical  extent  of  the  radar  beam.  The  effect  of  pulse  length  on  the 
filtering  of  the  gradient  is  negligible  since  the  150-m  range  resolution 
covers  only  an  8-m  altitude  change  at  3°  elevation  angle. 

Profile  3  (Figures  4.7  through  4.9)  again  illustrates  the  problem 
of  measuring  large  gradients  at  long  range. 

Profile  4  (Figures  4.10  through  4.12)  contains  large  low-level 
gradients  (as  might  be  expected  from  a  thunderstorm)  as  well  as  discontinuities 
at  mid- level  and  at  the  top  of  the  profile.  The  distance  or  height  interval 
perturbed  by  the  discontinuity  increases  with  range  or  height  because  of 
filtering  by  the  antenna  beamwidth. 

Profile  5  (Figures  4.13  through  4.15)  again  illustrates  the  ability 
of  the  radar  to  measure  wind  profiles,  except  near  discontinuities  or 
large  gradients. 

In  all  of  these  profiles,  the  results  for  a  6-degree  elevation  angle 
are  approximately  the  same  as  for  3  degrees.  The  major  difference  is 
that  for  a  6-deqree  elevation  the  minimum  altitude  is  doubled  (46  m) 
and  the  ranqe  at  maximum  altitude  is  halved,  so  that  large  gradients  are 
not  filtered  as  much.  If  there  are  significant  horizontal  gradients 
present  in  the  wind  field,  the  radar-measured  radial  velocity  profiles 
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would  be  different  for  various  elevation  angles.  Measurement  of  the 
profile  at  the  3-degree  and  6-degree  elevations  can  be  used  to  test  for 
important  horizontal  gradients. 

Model  6,  a  thunderstorm  downburst  model,  contains  horizontal 
gradients  of  the  horizontal  wind  that  show  the  problems  associated  with 
radar  location.  The  figures  are  also  In  groups  of  three  with  the  same 
display  as  used  above.  Results  are  shown  for  an  Input  profile  that  places 
the  location  of  the  downburst  near  the  touchdown  point. 

In  Figures  4.16  through  4.18,  the  radar  location  Is  at  the  touchdown 

point  with  the  downburst  offset  707  m  (X  ■  Y  =  500  m).  Figures  4.16  and 

c  c 

4.17  show  that  the  radar  can  measure  the  profile  when  the  radar  beam  Is 
along  the  3-degree  glide  slope;  however  when  the  antenna  elevation  angle 
Is  6-degrees  (Figure  4.18),  the  radar  observes  a  significantly  different 
part  of  the  downburst  at  every  height,  and  the  measured  profile  deviates 
from  the  actual  profile  along  the  glide  slope. 

Figures  4.19  through  4.21  show  the  results  when  the  radar  is  displaced 
by  100  m  from  the  touchdown  point  in  both  x  and  y  (XR  3  VR  *  -100  m), 
while  the  downburst  is  centered  at  Xc  ■  Yc  *  +500  m.  The  3-degree  glide 
slope  profiles  (Figures  4.19  and  4.20)  differ  only  slightly  from  the  actual 
profile,  while  the  6-degree  elevation  angle  profile  deviates  significantly 
as  before.  Displacement  of  the  radar  to  XR  ■  YR  ■  +100  m  (profiles  not 
shown)  also  shows  a  small  deviation  of  the  radar-measured  profile  from  the 
actual  profile. 

However,  when  more  realistic  radar  offset  positions  are  used 
(X^  *  -2.5  km,  Yr  ■  0.1  km),  the  radar-measured  profiles  (Figures  4.22 
through  4.24)  do  not  even  resemble  the  Input  profiles.  At  every  height, 
the  radar  observation  is  made  several  kilometers  from  the  glide  slope  and 
the  result  Is  uancceptable  for  gradients  expected  In  thunderstorms.  For 


28 


offset  positions  this  large,  the  radar  would  have  to  be  scanned  in  azimuth 
or  elevation  or  both  so  that  the  radar  measurements  would  be  close  to  the 
aircraft  track.  When  the  radar  antenna  pointing  direction  Is  changed  so 
the  radar  pulse  volume  Intersects  the  glide  slope,  the  radar-measured 
velocities  approximate  the  headwind  encountered  by  the  aircraft  unless 
the  angle  between  the  radar  beam  and  the  aircraft  track  becomes  too 
large. 

When  the  horizontal  gradients  can  be  neglected  (as  in  profiles 
1  through  5),  both  headwind  and  crosswind  can  be  measured,  either  by 
rotating  the  antenna  through  a  complete  revolution  (the  VAD  technique) 
or  oy  measuring  the  profile  at  two  angles,  (Figures  4.25  and  4.26).  The 
wind  input  is  profile  2,  which  has  crosswind  as  well  as  headwind.  The 
antenna  elevation  anqle  is  3  degrees.  Figure  4.25  shows  the  measured 
profile  for  an  azimuth  angle  30  degrees  counterclockwise  from  the  glide 
path;  the  wind  direction  is  perpendicular  to  this  azimuth  anqle  for 
altitudes  greater  than  400  feet.  Figure  4.26  shows  the  profile  for  an 
azimuth  anqle  30  deqrees  clockwise  from  the  glide  path.  At  each  height, 
the  headwind  is  the  sum  of  the  two  radar  measurements  divided  by  2cos(30°), 
and  the  crosswind  is  the  difference  of  the  two  measurements  divided  by 
2sin(30°).  (Note  that  in  Figure  4.26  the  radar-measured  points  are  off- 
scale  and  not  properly  represented  between  250  and  450  m  altitude.) 


5.  CONCLUSIONS 


The  simulation  program  i  .  described  and  listed  in  the  Appendix. 

It  can  be  used  to  find  the  exiected  radar-measured  radial  velocity  profile 
for  any  antenna  position,  beamwidth,  pulse  length,  or  radar  location. 

The  wind  input  can  be  any  analytical  model.  At  present,  the  only  out¬ 
put  is  the  radial  velocity  profile,  but  the  program  could  be  expanded 
to  test  scanning  strategies  toat  would  combine  the  outputs  from  various 
antenna  pointing  directions  to  deduce  wind  profiles. 

The  simulation  results  show  that  for  an  antenna  beamwidth  of  1.5 
degrees  low-level  gradients  in  realistic  wind  profiles  will  be  properly 
measured  by  the  radar.  If  the  wind  field  does  not  contain  significant 
horizontal  gradients,  the  elevation  angle  is  not  critical.  Ihe  results 
Suggest  that  one  way  to  test  whether  horizontal  gradients  are  important 
is  to  measure  the  profile  at  various  elevation  angles.  At  an  elevation 
angle  equal  to  the  nlide  slope  (nominally  1  degrees),  stronn  oradients 
at  a  500-m  altitude  will  be  heavily  filtered  by  the  radar  pulse  volume, 
because  of  the  vertical  extent  of  the  antenna  beamwidth.  Therefore, 
this  stronq  filtering  will  remain  until  the  elevation  angle  is  very 
high,  and,  since  the  minimum  range  of  a  pulse  radar  is  at  least  several 
microseconds,  it  will  not  be  possible  to  resolve  strong  gradients  at 
altitudes  above  several  hundred  meters.  The  antenna  beamwidth  filters 
the  wind  profile,  but  does  not  siqni f icantly  filter  realistic  wind 
profiles  at  low  levels. 

The  effect  of  radar  offset  from  the  runway  presents  a  serious  problem 
for  wind  fields  with  significant  horizontal  gradients  coimonly  found  in 
thunderstorms.  It  is  evident  from  these  simulations  that  the  profile 
along  a  path  parallel  to  the  glide  path  is  not  acceptable  for  offsets 
of  several  kilometers.  If  there  is  a  precipitation  echo,  the  entire 
region  may  be  scanned  (rapidly)  and  a  single  Doppler  radar  may  be  able  to 


identify  hazardous  areas  without  actually  measuring  the  wind  vector.  For 
example,  a  contoured  display  of  the  radar-measured  radial  velocity  (obtained 
from  data  acquired  during  an  azimuth  sector  scan)  would  reveal  the  location 
of  the  downburst  in  model  6.  Figure  5.1  shows  the  radial  velocity 
profiles  measured  by  a  radar  scanning  a  downburst  at  several  azimuths. 

In  this  case  the  downburst  (model  6)  Is  located  2.5  km  from  the  runway 
(Xc  *  2.5  km,  Yc  *  0),  while  the  radar  is  located  at  XR  *  -2.5  km,  -0.1  km. 
The  elevation  angle  is  only  1-degree,  so  the  radar  can  detect  the  strong 
outflow  that  occurs  in  the  lowest  200  m. 
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APPENDIX:  COMPUTER  PROGRAM 


1.  Program  Structure 

The  computer  program  consists  of  the  modules  shown  in  the  structural 
diagram  (Figure  A.l).  This  chart  illustrates  the  heirarchy  of  the  modules 
and  not  the  sequence  of  execution  (Mvers,  1 975 ) “* .  The  connecting  lines 
indicate  that  the  module  below  is  called  by  the  one  above.  The  parameters 
passed  between  the  modules  are  listed  in  the  input-output  table  (Table  A.l). 
The  interface  numbers  refer  to  Finure  A.l.  The  parameters  listed  under 
INPUT  are  passed  from  the  cal  linn  module  to  the  one  called,  and  those 
listed  under  OUTPUT  are  returned  by  the  called  module.  The  parameter 
that  determines  which  wind  model  is  used,  MDN,  and  the  graphics  parameters 
are  the  only  ones  not  included  in  Table  A.l.  MDN  is  passed  by  means  of  the 
common  block,  MODEL. 

The  system  dependent  Graphics  routines  (namely  AUTOPIT, DDINIT,  and 
DDEND)  referenced  in  the  pronram  listinn,  and  their  associated  common 
blocks  (LABPL  and  SIZEPl)  have  been  omitted  from  this  discussion. 

2.  Description  of  the  Proqram  Modules 

GUS1M  obtains  the  wind  model  number  and  controls  the  computation 
of  the  wind  prof i les. 

GPPROF  computes  the  actual  radial  velocity  alona  the  glide  path. 

RMPROF  computes  the  radial  velocity  profile  measured  bv  the  radar. 

RPLOT  plots  the  actual  and  measured  radial  velocity  profiles  versus 
ranqe. 

ZPLOT  plots  the  radial  velocity  profiles  versus  heiaht. 

READGP  obtains  the  alide  path  parameters,  namely,  the  elevation  anale, 
the  azimuth,  the  maximum  ranqe,  and  the  ranoe  increment  at 
which  the  velocity  is  computed. 


"  G.  J.  Myers,  Reliable  Software  Thrcuah  Corr,K>8rte  Deeian,  (New  York: 
Petrocelli  Charter,  1975.) 
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READRM 


RNIflT 


EUNT 


AZINT 

RADVEL 


WINDN 


obtains  the  following  radar  parameters:  elevation  angle, 
azimuth,  maximum  runqe,  ranqe  increment  (for  integration  net), 
radar  location,  pu  se  length,  ranoe  delay,  and  beam idth. 
computes  the  radar  measured  mean  radial  velocity  within  a 
pulse  volume  by  integrating  in  ranoe  the  mean  radial  velocity 
In  the  beam  cross  section. 

computes  the  mean  radial  velocitv  in  a  beam  cross  section  by 
integrating  with  respect  to  elevation  the  mean  radial  velocity 
in  the  azimuth  beamwldth. 

computes  the  mean  radial  in  the  azimuth  beamwidth. 
converts  the  spherical  coordinates  to  Cartesian  and  computes 
the  radial  component  of  u,  v,  and  w. 

computes  the  wind  components  u,  v,  and  w  at  the  point  (x,y,z) 
for  wind  model  n. 


The  radar  antenna  illumination  function,  w  -  wawAZ,  where 

WfL  »  exp('^-'L--^--)  , 

2A2 


WA7  *  • 

AZ  2aZ 

was  factored,  so  that  the  integral  could  be  evaluated  one  step  at  a  time, 
as  follows: 

v  .?«“« 

, 

vei 


V  *  JL  y  V 

*M  Nl  VEL 


where  N  Is  the  number  of  range  Increments  in  a  pulse.  This  was  done  to 
facilitate  proqraminq,  but  had  the  added  benefit  of  eliminating  the 
need  for  large  arrav  storage. 
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Figure  A.l  Structure  diagram  for  the  simulation  program. 


3.  Data  Input-Output 


This  program  has  been  n  n  In  conjunction  with  an  interactive 
program  that  interrogates  the  user  to  request  parameters.  The  parameters 
are  stored  in  a  file  named  TAPE23.  This  proqram  running  in  batch  mode 
accesses  the  file  to  obtain  the  parameters.  Many  of  the  input/output 
statements  are  list  directed  ( i . e . .  the  data  type  determines  the  format). 
To  improve  the  portability  of  the  proqram,  the  references  to  unit  23 
could  be  changed  to  a  standard  input  unit  and  the  list  directed  statement 
changed  to  formatted  statements. 

Except  for  listing  the  Input  parameters,  the  output  is  entirely 
graphical.  The  actual  glide  path  profile  is  plotted  as  a  solid  curve 
and  the  radar-measured  profile  is  plotted  as  one-point  per  pulse  volume. 
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4.  Listing 


INTER¬ 

FACE 

INPUT 

OUTPUT 

1 

GPW  Glide  path  radial  velocity 

GPG  Distance  alona  glide  path 

GPZ  Altitude  alonq  glide  path 

NGP  Number  of  points 

2 

RflW  Radar  measured  radial  velocity 

RRG  Radar  ranqe 

RMZ  Altitude  alonq  beam 

NRM  Number  of  points 

GPW.GRG.NPG.RMW.RRG.NRM 

4 

GPU.GPZ.NGP.RMW.RMZ.NRM 

5 

GPH  Glide  path  azimuth 

GPA  Glide  path  elevation  angle 

GPR  Maximum  distance  along  glide  path 
DRG  Distance  between  points 

6 

R  Same  as  GRG 

AZ  GPH  converted  to  radians 

EL  GPA  In  radians 

Radial  velocity  along  glide  path 

7 

RAZ  Radar  azimuth 

REL  Radar  elevation  angle 

RNG  Maximum  radar  ranqe 

DRR  Ranqe  increment 

XR.YR  Radar  coordinates 

RPl  Radar  pulse  length 

RRD  Radar  ranqe  delay 

RBW  Radar  beanwidth 

8 

RRG 

DRM  Pulse  length  in  km. 

ORR 

AZ  RAZ  in  radians 

EL  REL  In  radians 

XR.YR 

BW  RBW  in  radians 

Mean  radial  velocity  in  pulse  volume 

9 

R  Same  as  RRG 

DRR.AZ.EL.XR.YR.BW 

Mean  radial  velocity  In  beam 
cross  section 

10 

r.orr.az.el.xr.yr.bw 

Mean  radial  velocity  with  respect 
to  azimuth  beamwldth 

11 

R.AZ.EL  ,YR,YR 

Radial  velocity  at  a  point 

12 

X ,Y,Z  Coordinates  of  point 

U.V.W  Wind  components 

Table  A. 1 
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PROGRAM  GUS  ! M (  l  NPU1  OUTPUT  ,  TAPEaM  ) 

COMMON /MODEL /MDN 

DIMENSION  GPW (  1001  )  ,GR(i(  1001  )  ,GPZ(  1001  )  , 
X  RMW(IOI)  ,RRG(10U,RHZ(101) 

CALL  DDI  N  I  T  ( <> ,  <?0H  SWEEZY  X6«6<>  RB(  ) 

REWIND  <>3 

CONTINUE 

NRM  =  100 

NGP  a  1000 

RE  AD  ( «?  I ,  •  )  MDN 

PRINT*,"  MODEL  NUMBER  ",MDN 
IE  (MDN  .LT.  1  )G0  TO  ✓ 

IE ( MDN .CT.6)G0  TO  ^ 

CALL  GPPROE (GPW.GRG.GPZ  ,NGP) 

CALL  RMPROE ( RMW , RRG , RM7 .NRM) 

CALL  R PLOT (GPW,GRG,NGP,RMW,RRJ,NRM) 

CALL  ZPLOT (GPW.GPZ ,NGP, RMW.RMZ ,NRM) 

GO  TO  1 
CONTINUE 
CALL  DDFND 
END 


SUBROUTINE  GPPROE (GPW ,GRG ,0P7 ,NGP) 

ROUTINE  TO  COMPUTE  RADIAL  VELOCITY  PROFILE  ALONG  GLIDE  PATH 
GPW  IS  RADIAL  VELOCITY 

GRG  IS  RANGF 

GPZ  IS  ALTITUDE 

NPG  IS  NUMBER  OE  POINTS 

DIMENSION  GPW( 1) ,GRG( 1  )  ,  »FZ(  1  ) 

CALL  READGPfGPH.GPA.GPR,: RG) 

CONVERT  ANGLES  T  RADIAN 
AZ  *  0.017«S33*GPH 
FL  *  0.0174SM*GPA 
NG  a  GPR/DRG  ♦  1 
NGP  *  MINO(Nf,P  .  NG) 

R  *  0 

DO  5  K  i  1 , NGP 

R  a  R  ♦  DRG 

GRG(K)  a  R 
GPZ(K)  S  R  *S  I N  (EL ) 

GPW(K)  a  P  A  DVE  L ( R , AZ , EL, 0.0, 0.0) 

CONTINUE 

RETURN 

END 


o  o  o  o  o 


SUBROUTINF  RMPROF ( RMW , RRG , RMZ , N RM ) 

ROUTINF  TO  COMPUTE  RADAR  Mr  ASURFD  RADIAL  VFLOCITY  PROFILF 
RMW  IS  RADIAL  VFLOCITY 
RRG  IS  RANGE 
RMZ  IS  ALTITUDE' 

NRM  IS  NliMBF R  OF  POINTS 
DIMFNSION  RMW( 1 ) , RRG ( 1  )  , RMZ ( 1  ) 

CALL  READRM(RAZ ,RFL, RNG.DRR.XR .YR.RPL, RRD.RBW) 

C  CONVERT  ANGLES  TO  RADIANS 
AZ  s  0.017«533»RAZ 
EL  =  0.0174533»RrL 
BW  =  0.017«S33*RBW 
C  CONVERT  RANGE  TIMFS  TO  KM. 

DRM  s  0. 15*RPL 
R RG ( 1  )  s  0. 15#RRD  ♦  0.5«DRM 
RMX  s  RNG  -  O.S»DRM 
K  s  0 

1  CONTINUE 

K  x  K  ♦  1 

RMZ(K)  r  RRG(K  )*SIN (ELI 

RMW(K)  x  RNINT(RRG(K),DRM,DRR,A7,EL,XR,YR,BW) 

RRG(K*1 )  =  RRG(K)  ♦  DRM 
IF ( RRG(K ) .LT. RMX  .AND.  K.LT.NRM1G0  TO  1 
NRM  s  K 
RETURN 
END 


SUBROUTINF  RPLOT (GPW ,GRG , NGP , RHW , RRG, NRM) 

C  ROUTINE  TO  PLOT  ACTUAL  GLIDF  PATH  PROFILF  AND  RADAR  MFASURED 
C  PROFILF  VERSUS  RANGE. 

COMMON /LABPL/LABX(8),LABY(8), ID(8) 

COMMON /SI ZF PL /XLO, XUP, YLO, YUP, KDX, KDY .KSX.KSY 
DIMFNSION  GPW( 1 ),SCL(7),KDV(7) 

DATA  SCL/O.  ,5.  ,  m.  ,  15.  ,iO.  ,  30.  ,40.  / 

DATA  KDV/0,Ut3,6,8,8.8/ 

LABX(  A  )  x  mHPANGF(KM) 

LABY ( U )  s  1  OH WIND  SPFFD 

XLO  x  0.0 

XUP  x  10.0 

KDX  x  10 

WMX  x  GPW( 1  ) 

DO  1  Js^.NGP 

WMX  x  AMA X 1 ( WMX , ABS (GPW ( J  ) ) ) 

1  CONTINUE 
DO  2  Jz2,7 

IF  ( WMX .  LT.  SCL  ( J  - 1  )  )G0  TO 
YUP  x  SCL(J) 

YLO  «  -YUP 

KDY  x  KDV(J) 

2  CONTINUE 

KSX  x  KSY  x  1 

CALL  AUTOPLT(GRG,GPW,NGP  ,«i ,8, 1) 

CALL  AUTOPLT (  RRG ,  RMW  , NRM , «? ,  -58 , 0) 

RETURN 

END 


SUBROUTINE  Z  PLOT (GPW  ,GPZ  ,  NGP ,  RMW  ,  RM7. ,  N  RM) 

C  ROUTINE  TO  PLOT  ACTUAL  GULF  PATH  PROFILF  AND  RADAR  MEASURED 
C  PROFILF  VERSUS  HEIGHT 

COMMON /LABPL/LABX(R  )  ,  L ABY (ft) , ID(8) 

COMMON /SI ZE  PL /XLO , XUP , YLO , YUP ,KDX ,KDY ,KSX , KSY 
DIMFNSION  GPW( 1 ) ,SCL(7  )  ,KDV(7  ) 

DATA  SCL/O.  ,5.  , 10.  , 15. ,^0. , 30.  ,40. / 

DATA  KDV/0,4,8,6,R,ft,8/ 

LABY (4  )  =  10HHE IGHT'KM) 

LABX ( 4  )  =  1  OH WIND  SPFFD 
WMX  =  GPW ( 1  ) 

DO  1  JsZ, NGP 

WMX  r  AMAXHWMX,  ABS(GPW(J))) 

1  CONTINUF 

DO  2  J  =  2,7 

IFIWMX.LT. SCL(J-1))G0  TO  2 
XUP  =  SCL(J) 

XLO  ;  -XUP 

KDX  =  K  DV ( J ) 

i.  CONTINUF 

YLO  ;  0. 

KDY  s  5 
YUP  =  0.5 
KSX  =  KSY  =  1 

CALL  AUTOPLT(GPW,GPZ, NGP ,Z,8, 1) 

CALL  AUTOPLT  (RMW,RMZ,NRM,2,-58,0) 

RFTURN 

END 


SUBROUTINE  RF ADGP(GPH ,GPA , GPR , DRG ) 

C  GLIDE  PATH  PARAMFTFR  INPUT  ROUTINE 

COMMON /LABPL/LABX (8)  ,LABY(8) , ID(8) 

COMMON 'MODFL/MDN 

RE. AD (<>3,  *  )GPH  ,GPA  ,GPR  ,DRG 

PRINT*,"  GLIDF  PATH  HF.ADING-180  (AZIMUTH)  "  ,  GPH  ,  "DEG  .  " 
PRINT*,"  GLIDF  PATH  ANGLF  (FLFVATION  )  "  ,GPA  , "DFG . " 
PRINT*,"  MAX  RANGE  ",GPR,"KM." 

PRINT*,"  RANGF  INC  RFMFNT  " , DRG , "KM  .  " 

ENCODE ( 30, 1002 , ID)MDN ,GPH ,GPA 
1002  FORMAT  (  "MODEL •  I  <?•  GLIDF  AZ  ,  EL*<?F5 . 0 ) 

RFTURN 

END 


SUBROUTINE  READRM! RAZ , REL , RNG , DRR , XR , YR , RPL , RRD, RBW) 
C  RADAR  PARAMETER  INPUT  ROUTINE 

COMMON/LABPL/LABX ( 8 ) , LABY (8 ) , ID(8) 

READ(23,*)RA7.  ,REL,RNG,DRR,XR,YR,RPL,RRD,  RBW 
PRINT  1002 ,  RAZ 

1002  FORMAT! •  RADAR  AZIMUTH  ANGLE  *F4 . 0«DFG . • ) 

PRINT  1004, RFL 

1004  FORMAT ( •  RADAR  FLFVATION  ANGLE  •F3.0«DFG.«) 

PRINT  1006, RNG 

1006  FORMAT ( •  MAX  RADAR  RANGE  •F4.1«KM.») 

PRINT  1008, DRR 

1008  FORMAT ( •  RANGE  SAMPLING  INCREMENT  »F4.2#KM.«) 

PRINT  1 0 1 0 , XR  ,  Y R 

1010  FORMAT ( •  RADAR  DISPLACEMENT(X.Y)  FROM  TOUCHDOWN  (• 

X  F4.2»,«E4.2»)KM.») 

PRINT  1012, RPL 

1012  FORMAT ( •  RADAR  PULSE  LENGTH  »F 3 . 1 *MICROSFC . • ) 

PRINT  1014, RRD 

1014  FORMAT ( •  RADAR  RANGE  DELAY  *F3. 1 •MICROSEC .• ) 

PRINT  1016, RBW 

1016  FORMAT! •  RADAR  BEAM  WIDTH  "E 3 . 1 •DEG . • ) 

ENCODE (2  3,2002, ID(4 ) )RAZ  ,  REL 
2002  FORMAT! •  RADAR  AZ,EL»2E5.0) 

RETURN 

END 


FUNCTION  RN!NT(RRG, DRM.DRR, AZ.EL.XR ,YR,BW) 

C  COMPUTES  MEAN  RADIAL  VELOCITY  IN  PULSF  VOLUME  BY  INTEGRATING, 

C  WITH  RESPFCT  TO  RANGE,  THE  MEAN  VELOCITY  WITH  RESPECT  TO  BF AM  WIDTH 
R  =  RRG  -  0 . S*DRM 
RMX  s  R  ♦  DRM 
SV  :  0. 

EN  s  0. 

1  CONTINUE 

VR  =  ELINTfR, DRR, AZ.EL.XR, YR.BW) 

SV  =  SV  ♦  VR 
EN  *  EN  ♦  1 . 

R  =  R  ♦  DRR 
IF ( R . LE . RMX )G0  TO  1 
RNINT  =  SV/EN 
RETURN 
END 
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F  UNCHON  EL  I  NT ( R , DR3 , AZ , EL , XR , YR , BW) 

C  COMPUTF3  MF  AN  RADIAL  VFLOCITY  ACROSS  BF AM  BY  INTEGRATING, 

C  WITH  RESPECT  TO  FLEVATION  ANGLE,  THE  MEAN  VELOCITY  WITH  RFSPECT  TO 
C  AZIMUTH 

DELTA  =  0.3»BW 

ANG  =  3.1 •DELTA 

DEL  =  DRR / I R*COS ( EL  ) ) 

VL  =  AZINT(R,DRR,AZ.FL,XR,YR,BW) 

DL  :  0. 

bk  =  k .  •DELTA*DFLT  A 
SW  =  I. 

1  CONTINUE 

DL  :  DL  ♦  DFL 
ELP  =  FL  ♦  DL 
ELM  :  FL  -  DL 

VP  =  AZINT(R,DRR,AZ,FLP,XR,YR,BW) 

VM  =  AZINT(R,DRR,AZ,FLM,XR,YR,BW) 

WT  =  FXP(-DL»DL/D«>) 

SW  =  SW  ♦  k . *WT 
VL  s  VL  ♦  WT«(VP*VM) 

I F ( DL . LT . ANG)GO  TO  1 
ELINT  =  VL/SW 
RETURN 
END 


FUNCTION  AZINT(R,DRR, AZ.FL.XR ,YR,BW) 

C  COMPUTFS  MF  AN  RADIAL  VELOCITY  WITH  RESPECT  TO  AZIMUTH 
DELTA  =  0. 3«BW 
ANG  =  3H  •DELTA 
DAZ  =  DRR/R 

VL  =  R  ADVEL ( R ,AZ,FL,XR,YR) 

DA  =  0. 

bk  :  k.  •DF.LTA*DELTA 
SW  »  1  . 

1  CONTINUE 

DA  s  DA  ♦  DAZ 
AZP  =  AZ  ♦  DA 
AZM  :  AZ  -  DA 

VP  =  R ADVEL ( R,AZP,FL,XR,YR) 

VM  s  RADVFL(R,AZM,EL,XR,YR) 

WT  s  EXP  ( -DA*DA/D«?) 

SW  =  SW  ♦  k.* WT 
VL  s  VL  ♦  WT* ( V P^VM ) 

IF ( DA . LT . ANG 1G0  TO  1 
AZINT  s  VL/SW 
RETURN 
END 
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FUNCTION  RADVFL ( R  ,  AZ  ,F.L,XR  ,  YR  ) 

COHPUTFS  RADIAL  VELOCITY  COMPONENT  OF  U,V,W  FOR  A  GIVEN  WIND  MODEL 
COMMON /MODEL/MDN 
SA  s  SIN(AZ) 

CA  s  COS(AZ) 

SF.  =  SIN  (EL) 

CE  =  COS (EL) 

X  =  R*CE*SA  ♦  XR 
Y  s  R#CEBCA  ♦  YR 
Z  =  AMAX1 (O.O, RBSE  ) 

IF ( MDN . EO. 1 )CALL  WIND1 (X , Y , Z ,U , V , W) 

IF  (MDN  .EQ.  i)CALL  W IND.2  ( X  ,  Y  ,  Z  ,  U  ,  V  ,  W) 

IF ( MDN . EQ. 3 )CALL  WIND? (X , Y , Z , U , V , W) 

I F ( MDN .  EQ .  4 )C ALL  W IND4 ( X , Y , Z , U , V , W) 

IF ( MDN . EO . 5 )C ALL  WIND5 (X , Y , Z ,U , V , W) 

IF(MDN.EQ.6)CALL  WIND6(X, Y.Z.U.V.W) 
rtADVFL  =  ( U  BS  A  ♦  VBCA)BCE  ♦  W*SE 
RETURN 
END 


SUBROUTINE  WIND  1 ( X , Y , Z , U , V , W) 

C  THIS  VERSION  GENERATES  U  GIVEN  Z 

C  THIS  PROFILE  IS  FOR  NO  VARIATION  OF  WIND  DIRECTION  WITH  ALTITUDE 
C  .00045. LE. Z. LE. .45KM.  ( 1 . 5 . LE . Z . LE . 1 500F I . ) 

C  U<>0  IS  WIND  AT  «?OFT .  ALTITUDE 
C  NOMINAL  VALUE  s  5  M/S 

C  FOR  REFERENCE,  CONSTANTS  DEFINED  AND  FORMULA  GIVEN  IN  FEET 
C  IF( Z . LT . 1.5)  GO  TO  600 
C  DL0G1 3=LN  13.33 
C  ZLOG 1  5=LN  (Z/1.5) 

C  UsU<?0*ZL0G  1 5/DLOG  1  3 
U006s 1 1 . 5 
UsO 
WsO 
VsO 

IF(Z.LT. .00045)  GO  TO  600 
DLOG 1  3  =  <? .  5900 171  34 
ZL0G15=AL0G(Z/. 00045) 

UsU 006* ZLOG  1 5/DLOG  1  3 
RETURN 
600  UsO 

RETURN 

END 
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SUBROUTINE  WIND  <MX,Y,Z,U,V,W) 

C  THIS  VERSION  GENERATES  U,  V  EOR  GIVEN  Z 
C  WHERE  0.LE.Z.LE.0.45KM. 

C  NOMINAL  VALUE  OF  U006  IS  0  M/S 
C  THETO  IS  SURFACE  WIND  (.006  KM.)  DIRECTION 
THETO* 1 80 . 

UsO 
V  =  0 
WsO 

IF(Z.LE.O.O)RETURN 
C  FOR  WIND  SPEED  WE  HAVE  , 

C  U006  IS  WIND  AT  .006  KM. 

U006=3.5 

UU  =U  006  •(Z/.O06)**0.*<3 
C  FOR  WIND  DIRECTION  WE  HAVE, 

TF((.12.LE.Z).AND.(Z.LE..«5>)  CO  TO  14 
IF(<  .03.LE.Z)  AND.  (Z.LF.  .  1«?) )  GO  TO  15 
IF( (O.LE.Z) . AND. (Z.LE. .03))  GO  TO  16 
C  MODEL  FAILS  FOR  THIS  7. 

RETURN 

14  THFT=  1  <?0*THFT0 
GO  TO  17 

15  THET=1<?0.»(Z-.03)/.09  ♦  THETO 

GO  TO  17 

16  THF.TiTHFTO 

17  U*UU»COS(THET/57.3> 

VsUU»SIN (THFT/57. 3) 

RETURN 

END 
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SUBROUTINE  WIND  3  (  X  ,  Y  ,  Z ,  U  ,  V  ,  W) 

C  THIS  PROFILE  IS  THUILLIER/LAPPE  NUMBER  5 
C  FOR  NIGHTIME  STABLE  LAYER  BOUNDARY 

C  EXPECTS  DISTANCE  IN  KILOMETERS  GIVES  U  IN  METER/SECOND 
C  U009  IS  VELOCITY  AT  .0091  KM. 

DIMENSION  UZ( 12) ,Z0( 12) 

DATA  (UZ=1 . 389,2. 109,3.550,4.270,5.866,5.403,4.734,4.219, 
A4. 065, 4. 013, 4. 013) 

DATA  (Z0s.0091 , . 0228 , . 0457 , . 09 1 4 , . 1 37 1 , . 1 828 ,. 2286 ,. 2743 , 
A. 3657, .3962, .4572) 

C  FIND  HEIGHT  INTERVAL  I 
U009=UZ ( 1 ) 

U=V=W=0. 

IF ( Z . LE . 0 . 0 ) RETURN 
DO  80  J  =  1  ,  12 
K  =  J 

IF(Z.LT.Z0(J))  GO  TO  85 
80  CONTINUE 

C  MODEL  FAILS  FOR  THIS  Z. 

RETURN 
85  I  =K 

IF ( I .GT. 1 )  GO  TO  3 

2  UsU009*(Z/.0091 )»«.45 
RETURN 

C  DO  LINEAR  INTERPOLATION  FOR  1=2,12 

3  UZIsUZ(I) 

UZ2=UZ( 1-1 ) 

Z01=Z0(I) 

Z02 sZO ( I  -  1 ) 

U=UZ2-(Z02-Z)«(UZ2-UZ1 )/(Z02-Z01 ) 

RETURN 

END 


4.116, 

.3200, 


4S 
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SUBROUTINE  WIND  4 ( X , Y , Z , U , V , W) 

C  THIS  PROF  ILF  IS  FOR  THUNUE  HSTORM  COLD  AIR  OUTFLOW 
C  EXPECTS  DISTANCE  IN  K I LOMETFRS  GIVFS  U,V,W  IN  MFTFR/SF COND 
C  THIS  PROFILF  GIVFS  W,  GIVKN  Z 

C  FOR  W,  THFRE  ARE  4  CASES;  MAJOR  AND  MINOR  UPDRAFTS  AND  DOWNDRAFTS. 

C  CONSTANTS  USFD  IN  THIS  MODEL  ARE;  ZD,ZR,A,L,P1,P2,Q0,Q1,Q2. 

C  VALUES  FOR  ZD.ZR.L  ARE  IN  KM., A  IS  IN  MFTER/SFCOND 
DIMENSION  VZ (7 ) ,UZ ( 7)  ,  Z0(7  ) 

DATA  <UZ=-5. 1329,-5. 1328,-2.5664,7.6992, 15.3984,3.5929,3.5929) 

DATA  (VZ* 1.5 398, 1.5398,4. 1062,4. 1062,5. 1328,6.6726,5. 1328) 

DATA  (ZOsO, .0061 , .0609, .0900, .1350, .2250, .4572) 

DATA  ZD.ZR.XL.A.PI ,P2, 00,01, Q2/. 2347, . 1371 , . 0396 , 4 . 2672 , 

41.21,0. 36,-0. 36,2.4165, 1.0769/ 

U=V=W=0 

IF(Z.LF.O. ORE  TURN 
C  FIND  HEIGHT  INTERVAL  I 
DO  80  Jsl ,7 
K  :  J 

IF ( Z . LT . Z0( J ) )  GO  TO  85 
80  CONTINUE 

C  MODEL  FAILS  FOR  THIS  Z 
RETURN 
85  IsK 

I F (  I  . GT . 1 )  GO  TO  3 

2  Us-5.1328 
V  =  1 . 5  398 
GO  TO  49 

C  DO  LINEAR  INTERPOLATION  FOR  1*2,12 

3  UZ1*UZ(I) 

UZ2=UZ(  1-1  ) 

VZ  1*VZ(  I  ) 

VZ2  =  VZ  f I  -  1  ) 

Z01=Z0(I  ) 

Z02tZ0(I-1) 

U=UZ2-(Z02-Z)*(UZ2-UZ1)/(Z02-Z01  ) 

V=VZ2-(Z02-Z)*(VZ2-VZ1)/(Z02-Z01) 

49  CONTINUE 
X  X  =  Z  /  X  L 
XR=ZR/XL 

IF( (Z.GT. .234). OR. (Z.LT. .058))  GO  TO  55 
IF ( Z.LT.  .  1001  )  GO  TO  52 
IF( Z.LT. . 140)  GO  TO  51 

C  FOR  MAJOR  DOWNDRAFT  WE  HAVE  W  FOR  ZD . GE . Z . GT . Z R 

50  Ws-P1*A*SIN(3.1415927*(XX-XR)/Q1)  J 

RETURN 

C  FOR  MAJOR  UPDRAFT  WE  HAVE  W  FOR  Z R . GE  .  Z . GE . Z R-L 

51  WN= (A#(XR-XX)/Q0)*( (2*0 0^1 ) • ( ( X R-XX ) ••2- 1 ) - ( 3 •Q0»«2- 1 ) • ( X R-XX* 1 ))  j 
WD=((2»Q0*1 )»(Q0**2-1 )-(3*QO##2-1 ) • ( 0  0— 1 )) 

WsWN/WD 

RETURN 

C  FOR  MINOR  UPDRAFT  AND  MINOR  DOWNDRAFT  WE  HAVE  W  FOR  Z  BOUNDED  BY: 

C  MINOR  UPDRAFT,  Z  R- ( 1 4Q2 ) #L . GT . Z . GE . ZR- ( 1 42*02 ) *L 
C  MINOR  DOWNDRAFT,  Z R-L . GT . Z . GF . Z R- ( 1 4Q2 ) *L 

52  Ws-P2*A*SIN ( 3 • 1415927*(XR-1-XX)/Q2) 

RETURN 

C  Z  OUT  OF  RANGE 
55  WsO 

RETURN 


SUBROUTINE  WIND  5 ( X  ,  Y  ,  Z ,  U  ,  V  ,  W) 

C  THIS  PROFILE  IS  FOR  NIGHTIMf  STABLE  LAYER  BOUNDARY  (SEVERE) 

C  PROFILE  FOR  NO  VARIATION  OF  WIND  DIRECTION  WITH  ALTITUDE. 

C  EXPECTS  DISTANCE  IN  KILOMETERS  ,  GIVES  UZ  IN  METERS/SECOND. 

C  UZ  AND  ZO  ARE  TABLE  VALUES  FOR  LINEAR  INTERPOLATION 
C  U009  IS  VELOCITY  AT  .0091  KM. 

DIMENSION  UZ  (  1 «?  )  ,  ZO ( 1 ^ ) 

DATA  (UZ  =  3. 952, 5. 954, 10.060,12. 164, 16.630, 15.398, 13.  *»47, 12.626, 
AH  .754 , 1  1 .600, 1  1.446,1  1.446) 

DATA  (ZOs.0091 , .0228, .OH57, . 09 1 4 , . 1 37 1 , . 1 828 , .2286, .2743, .3200, 
A. 3657,  .3962, .4572) 

U009SUZ( 1 ) 

U=V=W=0. 

IF(Z.LE.O.O)RETURN 
C  FIND  HEIGHT  INTERVAL  I 
DO  90  J  =  1  ,  1  2 
K  s  J 

IF(Z.LT.Z0(J))  GO  TO  95 
90  CONTINUE 

C  MODEL  FAILS  FOR  THIS  Z. 

RETURN 
95  IsK 

IF ( I  .GT. 1 )  GO  TO  5 

4  UsU009*(Z/.0091 )  •• . 45 
R  F  T  U  R  N 

C  DO  LINEAR  INTERPOLATION  FOR  1=2,12 

5  UZ1=UZ(I) 

UZ2=UZ ( I -  1 ) 

Z01 =Z0( I ) 

Z02  =  Z0( I -  1 ) 

U=UZ2-(Z0<?-Z)#(UZ2-UZ  1  )/(Z02-Z01  ) 

RETURN 

END 


SUBROUTINE:  WIND  6(X,Y,Z,U,V,W) 

C  ROUTINE  TO  COMPUTE  A  DOWNRURST  WITH  CENTER  AT  (XC.YC), 
C  SCALE  HEIGHT  H,  AND  WIND  SCALE  WO. 

C  THE  FUNCTION.  WHICH  OBEYS  CONTINUITY  IS 
C  VR=WOB(A/H)BB.?BZBEXP(-ZBB,?/HBB,?)B(  1  -FXP  ( -  ( R /  A  )  )  /R 

C  W  =  -WO*(1-EXP(-(Z/H)»»«?:»FXP(-(R/ A  )••<?) 

C  LET  A  BE  .5  KM,  H  BF  . 1  KM,  WO  BF  6  M/S 
WO=6 . 

A=  .  5 
XC  =0 . 5 
YCsO. 5 
H= .  1 
DX=X-XC 
DY=Y-YC 

R<?sDXBB«?4DYBB«? 

IF(R^.EQ.O)  10, 20 
10  ERal. 

R=0 . 

GO  TO  30 

20  ER  =  EXP(-R<?/A#*«?) 

R  =  SQRT  (  R<? ) 

30  CONTINUE 

EZ=FXP(-(Z/H)BB<?) 

IF(R.EO.O)  40,50 
40  U=VsO 

GO  TO  60 

50  RF  ACT* ( 1 . -ER )/R 

VR=WOB(A/H)BB^BZBEZBRFACT 
U=VRBDX/R 
V  sVRBDY / R 
60  CONTINUE 

Ws-W0B(1 .-EZ)BER 

RETURN 

END 
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